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LEST OES PRCURES 


A Single Active Component With One Spare 


LZ 


I. INTRODUCTION 


In making a mathematical model for a real-life phenom- 
enon it is always necessary to make certain simplifying 
assumptions so as to render the mathematics tractable. One 
of the simplifying asSumptions that is often made is to 
assume that certain random variables are exponentially 
distributed. The reason for this is that the exponential 
distribution is relatively easy to wort sen. 

The property of the exponential life distribution which 
makes it easy to analyze is its (memoryless) lack of deteri- 
oration with time. By this we mean that if the life time of 
an item is exponentially distributed, then an item which has 
been in use for a certain amount of time is as good as a new 


item in regards to the amount of time remaining until it 


fa duleoe When the life time of an item is exponentially 
distributed, the failure rate function for the item is 
Const anes 


Under the assumption of constant component failure 
rates, it is possible to build a reliability shorthame 
[Ref. 1] for any system. The term system is used to desirupe 


a set of components organized to perform some mission. 


Any study on system reliability requires a description 
for the system's life anda derivation of the system's 
Survival function. The reliability shorthand gives a simple 
and easy way for describing a system's life, but Le = rs 
@ifficult to implement computationally since it involves 
considerable complexity in handling convolutions. Here, the 
term convolution refers to the summation of independent 
random variables (lives). 

This paper presents an analytical approach for obtaining 
a general equation for the survival function of any convolu- 
tion of independent and exponentially distributed random 
variables. 

Section 2 deals with convolutions in detail and gives a 
mathematical derivation for a general equation for the 
Survival function for any convolution of exponential random 
variables, using Laplace transforms and Theorem of ReSidues. 

Section 3 gives a mathematical derivation for an alter- 
native formula for the Survival function of any convolution 
of exponential random variables by computing the coeffi- 
cients of all the polynomials that accompany the exponential 
terms in the survival function equation. 

Appendix A contains a computer program written in 
Fortran for computing the reliability of any convolution of 
inependent and exponentially distributed random variables. 


This program uses the general equation in Section 2. 


Appendix B contains another computer program written in 
Fortran to compute the reliability of any convolution Wem 
exponential random variables, using the general formula in 


Section 3. 
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IIT. CONVOLUTIONS 


Bader “the assumption of constant component failure 
rates, Ehasmscetdonemeprescents ai analytical approach for 
finding a general equation for the survival function of any 
convolution of independent and exponentially distributed 


manadom variables. 


A. SURVIVAL FUNCTION BY INTEGRATION 


Let us start with a simple system; a standby system 


having one active component and ome cold Spare component. 


The life time of the system is 


T = Ty + iby) 3 
where 

By ee eNO 

ee EXP ( Awan). 
and 


Ty; T9 are independent. 


An active component A is to complete a mission of dura- 
tion t, while a spare component S replaces the active compo- 


nent when it fails. Miceli ec  auratitenmmort the active 


Wa 


~<——_— ME SSS SO eee 
haus 0 S E 


a 
S 
= WORKING COMPONENTS ———— ee 


Figure 2.1 A Single Active Component With One Spare 


component is T, and the life duration dor the spare compas 
nent, if if is used, is I, (see Fig 2.58 


The system survival function is given by 


Ps (2) = Gee =) 

t 

2 Bape) + J Frtes) tas) ds 
0 
X ’ y(t 
e it + if an : Aye ds, 
0 
where 
f+ (s) = - dF y (s)/ds 


re: 


\ e 


II 


By integration, the survival function fer seme system as 


he -Ayt M eet 
o- d, Mia A, 


Fry (Ce) = 


EZ 


Now, suppose a system consists of one active component 
with constant failure rate dy and two spares having constant 
failure rates No and ), respectively. The life time of the 


system is 


i= T, ae TS, 


where 
T, ~~) EXP(A,) . 
i EXP (A » ) 
} 
ee aur EXPN = et 
and 


Wa ho, hg BEES independent. 
ive Shorthand notation for the system is 
EXP(\,) + EXP()\>) + EXP(\z) 


The survival function of the sustem is given by 


Ng 


Fy (t) = Fre p(t) f Fay (es) tre Ble) ds, 
0 


where 
“ae (t) = ae rAt, AL shat 
ear , 
Ao- dy ie 
= —S 
Fy, (t-s) sAgtt ) 


1S 


and 


fT +1, (s) = -dF + ,(s)/ds 


LS | ~A2$_ eS) 


By integration, the survival function for the system is 


Fy (t) = Fret +%(t) 


h2d3 eM A, As € -Aet AiA2 € Ast 
Oe d0r \) (ATA2) (As 2) Oy aah) 


t 308 


To compute the survival function for a system of n 
lives, 


TS 1) seen oe 


we can proceed in a Similar manner. However the integrations 
will become increasingly complex, particularly if the tacit 
asSumption made so far that all component failure rates are 
different is abandoned. Therefore, we will use another math- 
ematical way to compute the survival function without direct 
integration. This will be done with Laplace transforms and 


the Theorem of Residues in the next SubSee@tmaons. 
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B. CLASSIFICATION OF CONVOLUTIONS 


Consider a system which has a life time 


peat = ceca Lea, 
See, 5---,l, are independent and exponentially distrib- 
uted random variables. This system can be classified as one 


of three possible convolution cases 


l. First Case 


The system has n dissimilar failure rates, 


ig e.0U6UsmUtC—C8 5 Tis 


, t A; sl st Pies.) 
2. Second Case 
The system has n identical failure rates, 
dr; = \, eo = eee”. «Ths 
3. Third Case 


The system has a combination of the first case and 
the second case, Some Similar and some dissimilar failure 
mates. 

We consider the third case to be of particular 
interest, since it better represents the general situation 
in computing the reliability of any convolution of exponen- 


EHally distributed random variables. 
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C. CHARACTERISTIC FUNCTIONS AND LAPLACE TRANSFORMS 


Consider a random variable T which has a distribution 
FunecE Lon seis) and a density function f+(t). The chai 
teristic function of I is defined, for any) real nuns mae 


(Re iemec loeb: 


N 


4 


+ co 


fut 
if e alo Se iE, 


- © 


f(u) 


iT] 


+ oo 


| _fut 1g Ges eis 


— 0 


If T is non negative (T 20) and exponentially distributed 


with a failure rate io then 


f —(u) | atti dt 
0 
2 fut Ei 
-{ e a e dt 
0 
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Nowe 1 L 


ai. . aul, 
where 
Ti ~ EXP(\;), ele 1% 
and T,,...,Ty are independent, then the characteristic 
function of the sum (convolution) T,+ ... + T, is 


Bu) - [fu + Th) | 


“6- 
=t 
e 
"O.) 
sp 
G 


NH 


pry, 


TT (x48 ni) 


L=! 


ime implies that the characteristic function for the sum 
Meee tT ly «1S just the product of the individual charac- 


merrstic functions. 


Ly 


When dealing with random variables which only assume non 
negative values, it iS more convenient to use Laplace trans- 
forms rather than characteristic functions. The Laplace 
transform of the random variable having distribution F+(t) 


is defined [Ref. 2] by 


Gas 


ec 


D +(s) 


= 57 
= dF + (t) 


oC 


co 


; Sy 
e tome re at) ei Gue 


0 


mlbvevete 35 ae Ck is the density function of T. This integra 


exists for a complex variable s=at+bi, where az0. 


Ii T,,.-..,;I, are independent random variables, then Ga. 
Laplace transform of the sum (convoliition) 2, ©... .+ i.e 
given by 


on 4 ae) | 
ef o° ir | 


P 7,8) 


P (s) 


1} 
CH 
= a | 
ue 
vn 
— 
SS 


xe 
Bas 
2 
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: yf P 7.8). 


ieers important to note that the Laplace transform uniquely 
determines the distribution. That is, in fact, there is a 
one to one correspondence between distribution functions and 
Laplace transforms. 

Now, we will consider the following question. Given a 
imme t ion P+l(s)s dg@es there exist ,a function f+(t) whose 
Laplace transform is Dr(s), aide tetemexi Sts, how can 1t 
be determined ? The answer to the first question is not 
always positive. In’ general, the function PD +(s) must 
satisfy some restrictions which will be discussed in the 
next subsection. If the answer to the first question is 
positive, then the inversion of the Laplace transform by a 
complex integral gives the answer to the second question; 
imi: 18, the density function for the random variable T is 


given by [Ref. 3] 


—| 


f(t) = L(+) 





X + 1eo 
] Sk 
= a IT} ul e : P +(s) ds, 
XH — 100 


i 


where the notation [br(s)) refers to the inversion of the 
Laplace transform and the path of the integral above is the 


straight line Re s=5%(& >0) parallel to the imaginary axis. 


D. LAPLACE INVERSION AND THE THEOREM OF RESIDUES 


As was stated in last subsection, if we have a Laplace 
transform Dr(s) Eousmal non-negative random variable T, then 
the Laplace inversion L [pris] gives the density function 
toe) = 

When dealing with a convolution of random variables, 


that is, T = T,+ ... + Ty, , where T,,...,1, are independemas 


the convolution density function is given by [Ref. 3] 





7 Senne Lor] 
A+ 100 
t 
: _ J . . Seetey am 
xX - 10 


where PD + (s) is the Laplace transform of the sum 
(convolution )= i eel 

To obtain the convolution density function, it is Gleam 
from the above integral that we have to evaluate a complex 
integral, but the Theorem of Residues reduces the evaluation 


of a complex integral along closed contours to passages to 


the limit and ditterenerartrons. 
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Suppose a function G(s) is holomorphic ina finite. 
domain C, the term holomorphic meaning that the function 1s 


analytic and single-valued. Then the integral 


} 
a $ o(s) ds 





evaluated along it's boundary, vanishes. If the function 
G(s) is not holomorphic in C, the integral may not be Zero. 
iemenere 1S in GC only one singularity at point a, the value 


of the integral is called the residue of G(s) at the 
Singular point a. 


The well known Residues Theorem [Ref. 3] is 


The integral of an analytic function along a 
closed contour is the sum of the residues at 
the singular points in the domain enclosed by 


mee contour, multiplied by 2 TT i. 
Now, let us define the analytic function G(s) to be 
st 
G(s) =e . Drs) , 


where Pr(s) is the Laplace transform of the sum (convolu- 


Eden ) (leat esas ji, then, daceording to the Theorem of 
Residues, the convolution density will be obtained as 
follows; 


ZA 


| 
WN 
~6. 
— 
A 
a 








Eee Gey) 
A+ 1 0 
| sf 
= S ds 
2mT 1 ic Pri ) 
X= | 00 
Xx + 10° 
= G(s) ds 
Cilla | is) 
A —1 co 
m 
oa (ami) Ris 
Coa ae 
L=! 
where m is number of singularities (poles), and Rj , 
i=l,...,m, is the residue of the analytete tumetion Gc — ome 


the singular point (pole) aj 


Therefore 


™m 
L=! 


that is, the convolution density function is just the sum of 
the residue of the analytic function G(s) at the singular 
pPOLDES# 

The result established by the last equation would not be 
of much value if it were not possible to compute the residue 
of the analytic function G(s) at the singularities directly 


and without evaluating any integral. Thus, let uS now 


ae 


discuss how the residue at a Singular point (pole) a of the 
holomorphic function G(s), can be evaluated. 

First the order of the singular point (pole) must be 
established by computing the limit, for s tending to a, of 
(s-a)G(s), Lae Gis jee oidmao on,  tntil we find a finite 
limit. The exponent which we find for (s-a) is the order of 
Beeeesinesular point a. [f the singular point(pole) at a is of 
order (multiplicity) n, the residue of the function G(s) at 


fs given by [Ref. 3] 


n-! 


» ea i J ((s-2) 6¢s))] . 


(n-!)! gs7! 





Mig-retore, the residue of G(s) at the singular point (pole) 


a;,is 


et 


“| a on d— ((s- a) '6¢5)) | 


S~> Qj 








where ages the order (mult ipinciey eure Che Sincular point 


(pole) ai, eal, a mM. 
By Substitution, the convolution density function is 
Biven by 
™ 


f(t) = » Ry 


Zi 


™ n.-' 


E plat fale) 


}= 


This implies that the convolution density can be obtained by 
limits and differentiations which are easier operations to 
deal with than the integrations. 

Now, we are ready to deal with the three different 
convolution cases which have been mentioned before. we will 
try first to obtain the convolution density for every case. 
Then, the survival functions will be easily obtained by 


integration. 


E. SURVIVAL FUNCTION OF A CONVOLUTION 


In this subsection, we will use the results obtained in 
Subsection D to find the survival function for allege. 
possible convolution cases which have been mentioned before. 
This will result ina general equation for the survival 
function of any convolution of independent and exponentially 


distributed random variables. 
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1. A System Having n Dissimilar Failure Rates 
Consider a system defined by the shorthand notation 


EXP(A,) + ss EXP(\ _ ) 


where 


eee eee i,j <1, ... ,n. 
The life time of the system is 
Tig 


feere T; is an exponentially distributed random variable 


with a failure rate Ai; teenie 1, , » In are 
independent. The Laplace transform of the convolution is 
given by 
~s(T +--+ Tn) 
Pr(s) = Ele 
Ti 
ai ( Ss 
ep (N Smt De; 
If we define the analytic function 
st 
G(s)=e . Pr(s), 

then, according to the Theorem of Residues, the convolution 


density is given by 


25 


7 oe ) Ri 


where nis the number of singularities (number of distinct 
failure rates), and R is the residue of the analytic func- 
tion C(s) at the Simeularepoimnt wpome) ay ee 1 = bye ck 5 ee 
It is clear from the Laplace transform function 
“ye(e) that there are n singularities, each of which has 
an order (multiplicity) ome; that is, nj = 8°iorei -=2aee 


rie Therefore, the residue of G(s) at the singular 


point (pole) _Xr; is given by 


n,—! 


R. as, [Eom 5. ((s-¢ Ap) 6s) 


= lim cs de) G(s) 


Say 


By substitution, the convolution density will become 


f(t) 


alias 
a 


i 
— 
}- 
cS 
me 
ne 
i 
o 
Va) 
~. 
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The survival function of the system will be obtained 


meetollows ; 


F(t) = (een 


iT 
—, 
Hh 
a 

—_ 

c 
ie 

C 


meer is the well known formula for the survival function 


for a system having n dissimilar failure rates. 


20), 


2. <A System Having n Identical Failure Rates 
Consider a system defined by the Shorthand notation 
EXxP(\) +... + EXP( A) (n times ) 
The life time of the system is 
Dt ee he 


where. Ts 5 sede) =e leneeremre ay 1S an exponentially distributed 
random variable with a constant failure rate r, and T, aaa 


Tn are independent. The Laplace transform of the convouwme 


tron 1s 


Pr(s) 


if 
Bae f= 
@ | 
ie 
Fh 
al 
fas 
tT 
~~” 
O: 
=o 


Then 


i 
o 
Ta 
a 
O- 
|| 
a 


G(s) 
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In this case the analytic function €(s) has®@one s#eular 
mommnt (pole) of order (multiplicity) n at the point = = 
According to the Theorem of Residues, the residue of 


the analytic function G(s) is given by 


N~ | 


2 tio [aon Ge((5+) z=) 





a a (n 1) Sori((5+) “aay 


n 


| | . Stee 
= fim [4 45 (e | 


Poe eA Ss 


Nant Xt 


~(n-1)! : | 


Therefore, the convolution density is given by 


f(t) R 
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which is the density function ofm@the Erlame (nm 


distributac or 
Now, the survival function 1s eivenauy 
Fr(t) = P(T>t) 


co 


fs (ae) echt 

it 

= [. Xe rn why 
a laa | 


parts yields 


du 


Integrating by 


: ja” 


-rt 
(2-1)! = 


Fr(t) = 


Ar 
ct 
\V 


which is the well known formula for the survival 


the Erlang (n r ) distripution. 


3. A Sytem Having Some Similar 


Railune Rate 
The shorthand notation for the system is 
ExP(),) + + EXP(), ) 
+EXP(X, ) + EXP(A,) 
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and Some 


n, times 
(n, times) 


(n> times) 


Ried) 


function “om 


Dissimi laa 


«oe © @® ® © ® © #® 8 © 


+EXP(Ap,) +... + EXP(Am) (nm times) 


miere m is the number of distinct failure rates, andn is 
Mmeemultiplicity of the i-th failure rate i aoe ene 


mie Laplace transform of the convolution is 


m = 


Pris) =] [ i) | 


a S + hj 


Then 
st 
G(s) = e . D+(s) 
mM 
ale 
(4 
> 5 | J(—Ai _) 
a Saye 
In this case, the system has m distinct failure rates, So 
that the analytic function G(s) faseem singewlarities. 


According to the Theorem of Residues the convolution density 


function is given by 


mM 
fi) 2;., 


j=l 
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where R is the residue of G(s) at the singular point (pole) 


LIN 1 [Sele =: ee 


Since nj; is the multiplicity for the failure wage 
vi = 1. 23eane it follows that the the residue of the 
analytic function G(s) at the singular point (pole) \ is 
given by 
1 nj 
Rj = lim ——— £—.((s+Xj) G(s) 
PG ta 
S—>—/, o 
} 
Ni; -! 
Arr ™m nj 
| i Sh ee ee 
= lim Fran = wel aai((s +Aj) é eo | 
_ = ; 
ees ay Ar ds Je! | 


Jen ae ne pus a 
re! eA) 


os [tt 4 

pie (1, 
rah “a n= I th -7; 

| Ga Mi hoes -d — seh) ] 
(n;-!)! s2-Ajb ds” jai 





Let D = d/ds. Then D =d / ds and 


mM oa 


6 (TIN) tim, fo eT Miser) | 





i (nN;-1)! S-—>- 


JZ 


beet us define the functron 


mM 


-n; 
Hi(s) = ual (s + Aj). 


=\ 
j#i 


Then 


nil Tt sf st nj -| 
D e H,(3)] 4 e [(o a H; (8) 


and 


.”) 


R; in 2)! im |e (on ty H6) 


' 
From the Binomial Theorem the term (D + t) can be repre- 


sented as 


gee ni -! n-wl-K K 
@ + "i as ( ; t D 


Therefore 

Yr" n; 

UTX; n,. st er. yi 
1 a — ies a a ~ \K t D CS) | 
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S aN, ee n;-\-K 
raat . “| (.)t 4, Ks) 


K= 


Let us define the term 








K 
O(Gigh)) = = Hi 6S) ] 
a ee 7s | 
: ds‘ un($ +i) Pe. 
Then 
nt ee fie = A 


mi (TIA) 2D oe é 


Since 


m 
er(t)= ) RR; , 

jst 
it follows that the convolution density is 


m n-! Clip K) n;-K~-| _\;t 


f(t) (TTD, SS \ Reine é See Ol 


j=| K=O 
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This equation is a general equation for the density function 


of any convolution of exponentially distributed random 


variables. 
Now, the 


Bol Lows ; 


Be(t) = p( T >t ) 


V1 
—, 
rh 
Cc 
ee” 
Ow, 
e 


(TSS 


j=l K=0 


.m n-l nj-k 


EEL 


j=l K=0 22! 


This is a general equation for 


Stew bya | 


GGini) 


kK! (n;-K~1)! 


C(1,k) : 


nj-K is 
i Sez 


fumettom cam be obtained as 


fo @] 
nj-K-! | mie 
| u e r du 


a ze 
(Ajit) : jt 


ce) (Ait) ty 
KI co Gai Pe 


the survival function of any 


convolution of exponentially distributed random variables. 
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In the next subsection, we will give three examples 
of the application of the equation. These examples represent 


all possiple convoltMezeonetorms - 


F. EXAMPLES ON FRE GENERAL EQUATION 


1. A Convolution of Dissimilar Failure Rates 


PR PRE mm em 


Using the general equation which was derived in the 
last subsection, this example illustrates how to compute the 


Survival function for a system defined by the shorthand 


notation 
ExP(),) + EXP(\,) + EXP()\;) 
where 
d, # ee. \, # Nena Mot Az 
Since all the failure rates are dissimilar, it follows thas 


the total number of Gestinct farlosesnaeoon 


The multip@icity of the first distumiet Parle cae d, 1s 
ee ee 
The multiplicity of the second distinct failure rate Ky ts 


oe = i 
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—- _—_———————————————— ee 








Gitemmultiplicity of the third distinct failure rate \ , us 
nz = il. 


The survival function of the system is given by 


ZS Ma eK af | 
(t) (TLE y = Crtty Ke oul ty =Mit 


F 7 
I £5 Pore eel kK 
={ od 
it pal as K! di! z-—1|)! 
= suml + sum2 + sum3, 
where suml is the sum at i=l, Sum2 is the sum at i=2 and 


Ssum3 is the sum at i=3. 
a eye twee mas only the value 0, and z has 
only the value 1. The value of the derivative term C(1,0) is 


given by 





C(1,0) = 


| 
4 
Oz 
S&S 
_|« 
oN 
iN 
a 
j 
a3 
c= 


i 
apa oner 
| 
yao 
iN 
ae 
<< 
ae 
es 





( hom dy )( Ns ,) 
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Therefore 


j—1 I-K 


ii 1, \ = a 
Se d, ho Meee a a 4 a zh 











=e 
Koceee | 
eGo ae g yt 
= dy MINS 0! d, ( 7 S 
. ee | 
| le: ( dy- AL) ( Ag- AL AY 
ows Pa 


OA) OA) 


Similarly, at i1=2, mng=i1, k has only thejvaltwems 


and z has only the value 1. The derivative term €(2)0 jae 


‘cai de) (Asm ey 


lode : Agt 
(A\- he) \3~ No) 


C( 250) — 
and 


sum2 = 


Also, at i=3, nz = 1, k has omiy the value 0) Sauce 


has only the value 1. The derivative term C(3,0) is given by 


| 
C(3,0) =————____——___—_—_— S) 


(Aim As Ag As) 
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and 


Iw NS a A3t 
(A\- Xz)CA2- As) 


Sins = 


By adding suml, sum2 and sum3, the survival function 


of the system is 


ett) = suml + stimZ + sum3 


_ Ret 
dp As € =e A, Ase rat Ai A2 é€ ‘ 


“ Oe Osa) (u- ds ad) Ga DeN want 


Pe fA Gonvolation of Identical Failure Rates 


en ee ee ee 


This example illustrates how to compute the survival 


function for a system defined by shorthand notation 


EXP( A) + EXP( A) + EXP( ) ). 


Lit wcomvolmwntem bas Only one distinct failure rate 
Poem miltiplicity 3. Therefore m=l and n=3. The survival 


function of the system can be obtained as follows; 


Ni-! n-K Z—-\ rt 
Fr(t) - wy y ys Sa (A; t) - ae. 
=, oe MS * (z-1)! 
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a& ~ ei) oe 
“RD K! oe (z-!)! 


Since m= 1, andm, = 3, it follows that) 2 -) anc 


0,1,2. When k = O the derivative term C(1,0) is given by 


0 


rae y 
6(1,0) =| 45° ba(e +h) |. n 





| 


-| ike hy) = ol 





Proey | =) 
=a 
Since | 
[ Les N) | = 


it follows that the k-th derivative 
ja" / ds “Ts+y) _ 2 i > 
re 
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Mec 1s, C(l wees 0,.for ke = 1,2. Now, “the survival fune- 


tion for the system will become 


3 Zz 
(1,0) (At) -At 
Frity=h ) me 8 


Z=t 
3 , 
_At 
= a os <a , t Seer 


Ve | 


which is the survival function equation of the Erlang( ee 


@rctribution. 
oe conve llubmenm Olgsimitar and Ditssimilar Failure 
Rates 


This example illustrates how to compute the survival 
mimection for a system having some similar and some dissim- 
ilar failure rates, using the general equation. 


Suppose we have a system defined by the shorthand 


motation 


EXP(),) * EXP(\,) + EXP(A, ) 


In this case, we have two distinct failure rates d, and 


X - The total number of the distinct failure rates is 


41 


The first failure rate \, has the multiplicity 
in ae 

The second failure rate No has the multiplicity 
ngsl 


The survival function equation of the convolution can be 


given as follows; 
oe ; : cci,K) (A 

x i, K : 

Fr(t) =(T hi), T1* a 


= suml + sum2 , 


where suml is the sum at i=l, and sum2 is the sum at i=2. 
At i=l, k takes the values O and 1. We will evaluate 


the two derivatives C(l@0) and C( 1p mecwirolows: 





Za 
0 TI el 
} f 


| Ts J] 


J#I 


hee. 
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2 
2 Tid” 
@r1,1) al Je! jg) (S + dj) ™ i, 
; =| 
O27.) 


Then 


suml = r x yy Spm K) Cat) oMt 


Qe 
Ks0 z= Ky MY (z-1)! 





2. CCUM -At cHraht-At cou) -At 

Pee | 

eat Adet dds a1 
hg~ Ayr Ae= Ay Ou Ap? 


At i=2, mp ewand k takes omily the value Oo.”  Thws, 
we have to evaluate only the derivative term C(2,0) as 


Pomeows ; 





Z, 
6) — 
a2,0) = E TACs + Ai) |. 
| 


4,3 


Then 


» One ae 
sum2 = hy hyn CCA, ) (A2t) 


0! ho 0! . 
2 | = hme 
\, As=ooos-—<C 
| “(dr Aa) AD 


(eo 


dy =e 
eco e 
Kee ) 


Now, the survival function for the system is given by 


suml + sum2 


F > (t) 
Sak 
ae : Ake ye" Ae t 
ho d, ho \, ro \, Ok) 


The three examples above illustrate how we can 
compute the reliability of any convolution of exponenial 
random variables. In practice, a computer program computes 
the reliability directly. A practitioner of the shorthand 
methodology can use the computer program in Appendix A to 
obtain a system's survival function by just inserting the 


failure rates of the system components. 
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IIL. SURVIVAL FUNCTION AND POLYNOMIAL COEFFICIENTS 

In section 2, we derived a general equation for the 
Survival function of any convolution of independent = and 
exponentially distributed random variables. This section 
gives an alternative formula for the survival function, 
using the idea of computing the coefficients of the polyno- 
mials that accompany the exponential terms in the survival 
function equation. This idea was explored by Sadan Gursel 
[Ref. 4], adding one distribution at a time. In this section 
we Will derive the same polynomial coefficients from the 


general equation developed in the previous section. 


Pee DETERMINATION OF THE POLYNOMIAL COEFFICIENTS 


Consider a system defined by the shorthand notation 


EXP(), ) eT Senet EXP( A, ) 1, times) 
9) Bae EXP(), ) (n, times) 
+EXP(A Ve As EXP(A _ ) (Ny times). 
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m 
This notation represents the convolution of \ n; exponen. 
=I 


tial random variables, where there are n identical exponen- 


tial random variables having the failure rate 


pod = 1, ...,m 


We expect the survival function of the convolution to 


have the form 


Sal ei ee A,\(t).e SA) 6 te eee , t > Om 
where 
2. il ae 
A,(t) = ayo tay tray, tt + ain \ E , 
7) nN2-! 
A,(t) = ag t a5, t+ agg t t..... Ao ry! E : 
+ 
a mo! 
Arm GE) = amo @ Sane. 2 ieee amn tt 
Them@metat vom fy (Ce i=l,...,m, represents a polynomial of 
(nj -l)st degree. Therefore, to obtain the survival funcouem 


of the convolution, we have to find all the coefficients 


5 El, « sie Sy j= 0 eee ee 


46 





For example, the convolution represented by the shorthand 


notation 


EXP(,) + EXP(\,) + EXP(\,) + EXP(X,) . 


where 


it ho 


has a survival function represented by 


BMGE)P oA, (e) . =A,e + Ao(t) che ‘dt 
where 
Z 
Ay(t) = aig + ay, t+ aygt ; 
and 
A,(t) = azo 


It is important to note that the number of exponential 
terms iS equal to the number of dissimilar failure rates, 
and each exponential term has a polynomial coefficients of 
degree one less than the number of identical random vari- 
ables having the corresponding failure rate. 

The problem of computing the coefficients will be solved 
if we derive a general formula for the constant coefficients 
of the polynomial coefficients in the Survival function 


equation. 
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Let us begin with the general equation which is 


m nj mM Nn) Nek . Z-! oe ic 
= Ci, k) CA;t) a - 
— CTA 2 ys sy eI! a 


7) ° 
{ ; 
Tne sera | ye: Can be written 
ixl i 


ee 
je a AO 
iri 


Substituting the above into the general equation yields 


Te 
, KE yp iS @-n! 


Z| 
m n: m n; Wel WK GG RCE) -r.¢ 
Fe) =) ACH ) d Fa 
i= = 


Sei Ke 


Rearranging the above equation yields 


_ z\it i “cli,k K+Z-1 Ze! 
F(t) = a TIN pi ae. cer ot 
=| ‘ai K=0 Z=| KY (- 
Let L = z- 1. Then the survival function will become 
m 
Ait L 
— eo e ne aor: ; 
laa i CHA; xy. 7 K! aN i 
i=l J#i K=0L50 
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Let us define 


G; (k,L) = 
Then, 


Nj-1 (nj-\)-K 


0 +) G, (x,t) = ¢ 


K-0 L=0 +C 


+G 


+G 
or, equivalently 


N;—) (nN; -1)-K 


3 YG; (ie, L) = G 


K=0 L=0 We 


+G 


+G 


a K+L 
GGians | 
KI Lt - 


= 


MOO) + 0 se es +G; (0,n, -1) 
0) ere +G; (1,n, -2) 
MOP eOD ee + Gy( 2), 3!) 


oo ¢ @ @ &@ &@® &@® &@®# &# 8 @ 


*2« «© © © & & 


ees e0) 


mo: 0') AJ ara 1G «(Ogee onle)) 


/(O,1)+ eee ee. 6 +G Cu) 


WOOL 2 ie Ghat = 32.) 


j (O,n,-1) 
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ni (n;-)-L 


) \ Ce Geil). 


L=0 Kz0 
Therefore 
j n=} (n-.})=L + Leal 
Foye “ADDY Y CGE CK) te 
i KI De 
sei nag 
mM 
_\.t 
- \ o7: Aste) 
i=l 
where 
n; =! : (n;-1)-L +> Leal 
CU, K) it 
Aca) = a. * ei e 
l i (TA) leg! 
nee Lt 
Now, it easily follows that 
a v (4j-1)- Lb by 
(0 SE hn 
Jz! K! ee, 
jee) 
represents the coefficient of t for the i-th polynomial @a 
the survival function equation, L - Oye. BS LS 
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where 





, ™m fai 
C(a2) -| = Ld(s 7) ' 


This formula is the general formula for the coefficients and 
it can be considered as an alternative formula for the 
Servival function of any convolution of exponential random 
variables. 

It is clear from the above formula and from the general 
Saeyelon in section 2 that the most difficult part in the 
reliability computations is the derivative term C(i,k), 
particularly when dealing with complex convolutions’ that 
have a long stream of failure rates with high corresponding 
multiplicities. We will give the the practical method for 


computing the term C(i,k) in Subsection C. 


B. EXAMPLE 
Consider a system defined by the shorthand notation 
EXP(),) + EXP(A,) + EXP(A5) + EXP(\,) 


Then, the survival function for the system is given by 


hot et 


= ert e L 
Bet) = Ay(t)se “'s AQ(®).e ee ha (Ce)) oe 


5) Il 


The first distimet wait eerate My has a multiplicityaam 
the second distinct failure rate No haspa multiplicity ie 
and the third distict failure rate dz has a multiplicity we 


Therefore 


in as ee 
1) =a 
no ol, 
eS 


The polynomial ie) will be obtained from 


n,~ | 
yy ig) = a pues tte ain) t 
SUD St eur Mme 1 ai l in the above equation yields 
A,(t) = ay - 
Since 
m nn: njpat—lb IE 
(Toy, ae 
nih © (cre ist ‘acme 
J #i ag 
1t follows tia 
5 nN; I-!l-0 K+0 
: C(!,K) 
| } 
aig = ( = Mi) K! 0! r 
j#! K=9 . 


By 





2 e K 
: hes + ¥ CCL AK) \, 


| 
K=0 St 


Za 
=o 3 Ci) 5) 


where 
2 nj 
C(1,0) = - Leen ee). d, 
; | 
Thus 


ay 
POG) O5-i)” 


Now, the polynomial ct) 1 


A, (t) = 410 


ho he 
eo A-A,)° 


Similarly, the polynomial Ag(t) 1S even (Dy 


Ay (t) = 290 
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2 
hy X3 
Z 
(Ay - A2)(As -X2) 
The polynomial Az(t) will be obtained from 
nz-! 
Az(t) = 25 e eooee e agn,-1 ¢ 
SUDSTLEULING Goo eclds 


Az(t) = az, * ag, t - 


In this case we have fo evaluate thewevowcocrtt ter ents 
a z9 and az; - substituting nz =2, m= 597 L= 0 ita 


general equation of the coefficients yields 


3 nj (ENN 8, K +0 
SES). 
azso=( TAG DD KI ol. 
J#3 K =0 e 
\ K 
Cis 2) 
= dy 2 >. “ Ms 
K=0 


=), eno) F dN Age(351) ie , 


where 
-n; 


[a lee) J. 


CiC32 0) 


0 = yee 


(Aim Ag (Ago As) 
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3 its, 
d 
con =[ Tse) Jeon, 


| | 
; | Ody As) ( 0; oe | 


aus 


A, Ag A, AzAz NK AoXds 


feeeeereutane t= 1, n,- 2, and m = 3 in the general equation 


of coefficients yields 








3 nN: Z-I-! K+ | 
J Cay K) 
asi = ( Any \j ) D Ki 1! As 
S93 K 30 
0 K+! 
es K) 
ak 1 | 
K=0 
= A, Az 6,9) do - 
Since 3 
4° ae 
(3,0) = | Ll (s+ Aj) | 
ds? ad ee 


J#3 


| 
( we XE hs) 


2 


it follows Vena 


 —————— 
(Ap ds)( ho ds) 
Finally, the survival function of the convolutionuees 
given by 
_ at ; ae 
Baie) = We. : ap(tyea ret . a3 
sag MT sag ehe™ + (agg rag) 802 


i» ene \, a Me 


( ho M( ee ) : ( AG ra) go 2) 


—=— 


i INS Io Ne = 
Oi ds) Xam As) CX Ag) OA Dg) 


As 7 
( Mr ds) Ne ds) ( Ni ere ds) 
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erie ner LGA METHOD FOR GOMPUTING THE DERIVATIVE TERM 


C(I,K) 


Smewot ene major ditficulties in computing the survival 
function either by the general equation, or by the general 


meemala of the coefficients, is the derivative term C(i,k), 


where 
K om =; 


C(i,k) = - TM(s +A;) ex,’ 


= 





KH cu 
1 fl 


A 
l 
} 
This term represents the k-th derivative for a product of 
m-l terms. The difficulty in comuting C(i,k) increases with 
k. Since k runs from zero up to ny gle, Pie: ie Lt ftoldows 
that the maximum value of k will not exceed the maximum 
Mrbeiplicity of the »*distinct failure rates minus one. For 
example, if we have the two distinct failure rates My and 
MS where \, hase aemtttapliciey mn, and > las a multiplicity 
No » Ny> No , then, the maximum value of k is hn ie 
In practice, due to the high accuracy (precision) 
required in the computations, we will limit the derivative 
term C(i,k) to values of k <9. This implies that the 
Meetaplicity of any of the distinct failure rates in any 


Somvolution must not exeed 10. 


a) Ti 


The following is one of the solutions which will be used 
in the computer program for computing the term C(i,k) (see 
Appendix. A, Appendix. B). 


Suppose we have a convolution which has m distinct 





failure rates with m corresponding multiplicities . We know 
that 
cask) = | T (seri) = N 
ds‘ Wei 7 l 
Suppose for example , we want to compute C(i,k) at i = l, 
then kiee0700 3 ae ies and 
-Ty 


C(1,k) 


ier ds*  Wftees) — 


| 
aa (Gang ed 


Our problem is how to change a derivative of a production 
terms to a derivative of a sum of terms, since the later is 
easier in computations. 


Let us define 


—_ TT (s+ hi) 


sey 


58 


he i-l we obtain 


wm”) - 7) 


| tg, (SN) 


Kh 
i 


; (S+ ho) (s+ 3) a (s+ OF a 


By taking the natural logarithim on both sides, we obtain 


N2 
liga f\= - [ints +A) Te eect) (Ss dad 


Now, by taking the first derivative on both sides 


obtain 
Ea £ | = 4 si 
| | ae 2 © + Eee 
| S+)> Bo 


mee equivalently 


S l= * os f 


Let us define the term 


5) 9) 


, we 


Then, 


j=! she 
y#i 
= eee. cis 
Then 
= } 
f, = az Ay f 
or, equivalently 
df, /ds = am AY ° f,. 


When i=l and k=l, the term C(i,k) is 


m = Ty 
ee | = TMs +5) hemes 
j# } 


But 
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se that 


e(l,i) = - A, . a Y J 
= a 


We can proceed in a similar manner to compute the other 


mmmataves ©(1,2),6(1,3),...,€(1,n -1) as follows; 


Hh \ 


I 
Q. 
~ 
QO. 
Y 


d/ds - A | . fi 


put 
a _ 
fy = - Ay. £), 
so that 
FZ 7 2 
iB = - Ay ier 15 Ay fy 
a 
= CD) | Ag 4] f 5 
where 
ye ae Ay 
A, = d/ds 
jz S + Mj 
JF! 


m 4 ae 


J 
a ne ( S + \; ) 
jFl 


Then, the derivative term C(1,2) is 


a Pye.) 
2 ae: 
Cry fos =A, ee] nh , 


and so on until we obtain the derivative term Saye ih) (see 


Appendix.A, Appendix.B). 
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APPENDIX A 


This appendix contains a computer program written in 
[ierran for the survival function of any convolution of 
independent and exponentially distributed random variables, 


uSing the general equation 


. m nel - eit 
ay - (TTR, yyy EV > ua 8 ee 


a 
j=) K=0 


A. FUNCTIONS AND SUBROUTINES USED IN THE PROGRAM 
1. FUNCTION LFAC 


This function computes any factorial required in the 


program. 
2. SUBROUTINE GAMMA 


This subprogram computes the incomplete gamma func- 


tion with the parameters( ); ,n;-k) and returns the result in 


the variable G, 
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where 


for 


(1) 


(2) 
(3) 


Re z-| 


- me) =rit 
“e Ss 


3. SUBROUTUNE eae rv 


This subprogram computes the k-th derivative C(i,k) 


k<9 and returns the result in the variable C. 


PROGRAM'S LIMITATIONS 


The multiplicity of any of the distinet failure rages 
must not exceed 10. 

The complijerepreersieon limit. 

There is no limitation on the number of distinct faite 


rates. 
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section Computer program 


$ JOB 


Cr) OO Oe Ga" Oa @ 


a. 


2 


20 
10 


as) 


ADA... . 


ieee mee 


ee eee @ 


ee es @# @# @ 


P-tEhvedmatinectatarlure rate. 
eae cis tiivet tallure rates. 


7eroduee of all the failure rates. 


Suave! fimetiom of the convolution. 
Time. 

Value of the derivative term C(i, ven) 
Incomplete gamma function. 
Mipieeenetiy Of thesia-th distinct 
failure rate. 


Piaaveor MUletertelties of the distinct 


failure rates. 
OEMS E LICE tatlure rates. 
Loop index. 


DOUBLE PRECISION LAMDA(20), LAMDAI,PRODCT,FBAR, 


ol .c,G 


INTEGER I,J,N,M,K,MULTI,MULTIP (20) 


ERENT , 


READ, 


PRINT. , 


READ, 


PRINT, 


READ, 


T 


M 


Srlemee enter the time of duration’ 
‘Please enter THE # OF dissimilar failure rates' 


SimcemnEde Values or dissimilar failure rates' 


(LAMDA(J) , J=1,M) 
PRINT, ‘Enter corresponding multiplicities 


READ, (MULTIP(J) ,J=1,M) 

PRINT, OUTPUT 

PRN , == 52> 

is. ike) al 

FORMAT(' TIME = ',F10.5 

PRINT, ' FAILURE RATE MEE EPEECILTY 
FBAR = 0.0DO 


PRODCT = 1.0D0 


DO 10 


ia 


iM 


WRITE(6,12)LAMDA(1) ,MULTIP(T) 
FORMAT (8X,F10.5,10X,1I3) 

= MULTIP(I) 
PRODCT=PRODCT*LAMDA (I) **MULTI 
DO 20 J=1,MULTI 


PUTT 


K = J 


= al 


LAMDAI = LAMDA(T) 

CALL GAMMA(MULTI ,K,T,LAMDAI,G) 

CALL DERIV(I,M,LAMDA,MULTIP,K,C) 
FBAR=FBAR+C*G/ (IFAC(K)*LAMDAI** (MULTI-K) ) 


CONTINUE 

CONTINUE 

FBAR = FBAR * PRODCT 
WRITE(6,13)T,FBAR 

Remi ( bx,'P(T >',F6.2,') = ',F10.5) 
STOP 
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30 


40 


30 


70 
60 


a 


END 

FUNCTION IFAC(N) 

JFAC=1 

IF(N.NE.0O) THEN DO 

DO 30 J=1,N 

JFAC=JFAC*J 

CONTINUE 

END IF 

IFAC=JFAC 

RETURN 

END 

SUBROUTINE GAMMA (MULTI ,K,T,LAMDAI,G) 
DOUBLE PRECISION T,G,LAMDAI 

INTEGER L,K,Z,MULTI 

G=.0D0 

L=MULTI-K 

DO 40 Z=1,L 

G=G+ (LAMDAI*T )** (Z-1)*DEXP(-LAMDAI“T) /IFAC(Z-1) 
CONTINUE 

RETURN 

END 

SUBROUTINE DERIV(1I,M,LAMDA,MULTIP,K,C) 
DOUBLE PRECISION A,F,C,B(9),LAMDA(10) 
INTEGER L,J,M,L,IK,MULTIP(10) 

A=.0DO 

F=1.D0 

DO 50 J=1,M 

IF(I.NE.J) THEN DO 

F=F* (LAMDA(J)-LAMDA(I) )***(-MULTIP(J)) 
A=A+(MULTIP (J) / (LAMDA(J)-LAMDA(I))) 
END IF 

CONTINUE 

IF(K.GT.1) THEN DO 

IK=K-1 

DO 60 L=1,IK 

B(L)=.0D0 

DO 70 J=1,M 

IF(I.NE.J) THEN DO 
B(L)=B(L)+(((-1)**L)*MULTIP(J)* 


*( (LAMDA (J )-LAMDA (I) )**(-L=-1) ) )*I1FPAG (i) 


END IF 
CONTINUE 
CONTINUE 
ENE 


IF(K.EQ.0)THEN DO 
C=F 

ELSE DO 

TP CK B@e 1) THEN DOG 
CS ( = 1)*A*F 

ELSE DO 
IF(K.EQ.2)THEN DO 
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C= ( (APP, )-B( 1) )*F 

EUSE DO 

IF(K.EQ.3)THEN DO 

Cee) (A343 A B(L)-B(2) )*F 
Hor, DO 

IF(K.EQ.4)THEN DO 
C= (A**4-6*A**2*B(1)+4*A*B(2)+3* 
*B(1)**2-B(3) )*F 
Eiok DO 

IF(K.EQ.5)THEN DO 
C=((-1)* as *5) + 10% (A**3 )*B(1)- 1LO* (A**2 )* 
“bh (2)-15*A*B(1)**2+5*A"“B(3)+10*B(1)*B(2)-B(4) )*F 
BoE BO 

IF(K.EQ.6)THEN DO 

C= (A**6-15*(A**4 )*B(1)+20* (A**3 ) * 
*B(2)-15*(A**2 )*B(3 )+6*A*B(4)-B(5)+45* 
7 (A**2 )* (B ( li ) ww” ) -60.D0*“A*“B ( 1)* 
ee > Beye 3 )-15*BC1 )**3+ 10“B(2 )**2 )*F 
Er DO 

IF(K.EQ.7)THEN DO 

C= ( ( So ) v (A**7 ) +2] KAKNORB ( iE) -35*AwNKGKB (eZ ) + 325 ww Awe 3 ve 
*b (3 )-21*A**2*~B(4)+7*A*B(5)-B(6)-1LO5S*A**3* 
*b(1)**2+210*A**2*B(1)*B(2)-105*A*B(1)*B(3)-105* 
*B(1)**2*B(2)-70*a*¥ 
*B (2 )**2+105*A*B(1)**3+21*B(1)*B(4)+35*B(2)* 

(3) )*F 
EEoeE DO 

IF(K.EQ.8)THEN DO 

ce (A**8 -J2QXAKNEXB ( 1 ) +56 *ARNOKB ( 2 ) -7O*A**K4GKB ( 3 ) +56%* 
WARK 3*B(4)-28*A**2*B(5 )+8*A*B(6)-B(7)+210* 

wKARWAKB ( iL ) w*ED-S6O0*A** 3h ( uf ) Ww 
*B(2)+420*A**2*B(1)*B(3)-168*A*B(1)*B(4)+28* 
*B(1)*B(5 )-420*Ax*2*B (1) **3+840*A*B(1)**2*B(2 )+280% 
ee eee 2108 (1) 2B 3) -2807A*B(2)*B(3)-280*B(1)* 
See ore (2) -B(4)435"B(3)~*2+105*B(1)*+4)*F 
Eboe DO 

IF(K.EQ.9)THEN DO 

C= (-A**9+36*A**7*B(1)-84*A**6*B (2 )+126*A**5*B(3)-126* 
WAKEAKB se ) oe SAEs le ( 5 ) =3 SE ( 6 ) + ON AB ( 7 ) <8 ( 8 ) _ 
ee 21230 G1) *B(2)-1260*%a**3*B (1)* 
(3) +7563 Peed I) B= 252A” “ie oo. We} (| bee) )) 
=C+1260%A%*3%B (1) 3 - 3780%A%*2*B(1)% 
“B(2)- Seer oi 2 2+ 1890*a* 
3). 13 c04ae19%n(2)*8(3)-378"R(4)* 
*B(1)**2-404%A*B (2)*B(4)+84*B(2)*B(5)-945%*a% 
(1) **4+2520%A*B (1) *B(2)**2+1260*B (1) **3"B(2)-315" a 
*b (3)**2-1260*%B(1)*B(2)*B(3)+126"B(3)*B(4)-280%B(2)%* 
C=C*F 

Ee Le 

END IF 

END IF 


6/7 


END IF 
END IF 
END IF 
END IF 
END IF 
END IF 
END IF 
RETURN 
END 

SENTRY 
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Section Examples subsection All failure rates are dissimilar 


Shorthand notation 


EXP(.3)+EXP(.32)+EXP(.4)+EXP(.6)+EXP(.62) 


Oueput 
TIME = 2.00000 
FAILURE RATE MU etree CLTy 

O-30000 i 
02322000 ip 
0.40000 iE 
0.60000 iL 
0.62000 i 

mor > «6©22.00) = O99 Sak] 
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subsection All failure rates are ident eae 


Shorthand notation 


EXP(2.1)"+ <.. =a (Ze) (10 times) 
Output 
TIME = 3220 
FAILURE RATE MULT Teeer ry. 
2.10000 10 
P(T > “35208 Ome5772 
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co ct homsemnemsciiniitar and Some dissimilar failure rates 


Shorthand notation 


+ EXP(6.5)+ ... + EXP(6.5) (5 times) 
+ EXP(4.5)+ ... + EXP(4.5) (5 times) 
fee (3.5)+ ... + EXP(3.5) (5 times) 
feexP(2.5)+ ... + EXP(2.5) (5 times) 
pee (l].5)+ ... + EXP(1.5) (5 times) 
Met. 50)+ ... + EXP(.50) (5 times) 
eae put 
TIME = 8.00000 
FALLURE RATE Mort PriicertrTy 

6.50000 D 

4.50000 5. 

3.50000 5 

2.50000 5 

50000 5 

0.50000 5 
more> 8.00) = 0.99894 
cms cms 
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APPENDIX B 


This appendix contains an alternative computer program 
written in Fortran for computing the survival function of 
any convolution of independent and exponentially distributed 
random variables, using the general equation for the coeffi- 
clients of the polynomials that accompany the exponential 


terms in the survival function equation 


a L 
Fry- Yo e€ ' ) ai t , teeg 
j= sew 
where 
m nj (nj-N-L Cli, k) K+L 
| iy K 
iL ' ( TT A; ) ». Mj 

ae KI OL 
jt eS, : : 


A. FUNCTIONS AND SUBROUTINES USED IN THE PROGRAM 


LL... Fumctton LPAG 


This function computes any factorial required in the 


program. 
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2. Subroutine COEFF 


itemsubprosuam computes the coeffient of t for the 
meen pOlynomidl! “in the “Survival function equation and 


returns the result in the variable COQ, 


where 
mM nj n,;-I-l Cli K) K + L 
ly 
CO = (TTA ) - ‘eo ul hd; 
jf K=0 


Ss. SUBROUTINE DERIV 


This subprogram computes the k-th derivative C(i,k) 


mepeee-9 and returns the result in the variable C. 


B. PROGRAM'S LIMITATIONS 


memetie multiplicity for any of the distinct failure rates 
must not exceed 10. 

(2) The computations must be within the computer precision 
limit. 


(ye No limitation on the number of distinct failure rates. 
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section Computer program 


S JOB 


ANAANAAgQAAADAAIAaAaG 


diel 


10 


1 


14 


S10, 
20 


1S 


LAMDAIL...-.1=th distinct fazlure rate. 

VN GIDYG, Ao Array of M dastinct “fai lGrepeaces, 

MULE Dae Multaplicity of the’ Peth dasst mace 
failure rate. 

MULTIP....Array of multwpaieatres oF Enemeatoee ar 
failure rates. 

PRODCT....Product of the distinct failure rates. 

BAR s crue Survival function Gf tne con oMituaeme 

Lice... 4 ae Time. 

OM Value of the derivative term C(i,k). 

COx.. 9 Coeffetiont of the polynomials. 

M conte eeanaeee # of distinct failure rates. 

Li, J Skee: 5 Loop index. 


DOUBLE PRECISION LAMDA(20), LAMDAT,PRODCT,FBAR, 


<P Gaee 


INTEGER 1L,J9N.M,K,MULTIQ MULTE) 

PRINT, “Please enter Ene timemor aura ron, 

READ joes 

PRINT, ‘Please Enter # of dissimilar failure rates' 
READ, M 

PRINT, ‘Enter the values of dissimilar failure rates' 
READ, (LAMDA(J) , J=1,M) 

PRINT, Enter corresponding multiplicities 
READ eC Meneses) d= lM) 

PRINT CULeUr 


t 


PRIENiE. 'sssorcs' 

WRITE(6,11)T 

FORMAT(' TIME = ‘,F10.5) 

PRINT, ' FAILURE RATE MU IPLEGL Ty 


FBAR = 0.0D0O 

DO 10 I = 1,M 

WRITE(6,12)LAMDA(1I) ,MULTIP(T) 
FORMAT(8X,F10.5,10X,I3) 

BO 20 1 = 1.M 

MULTI = MULTIP(I) 

WRITE(6,13)I 

FORMAT(1X, 'POLYNOMIALS COEFFICIENT FOR LAMDA',I2,' is :') 
DO 30 LL = 1,MULTI 

ye eee 

CALL COEFF(I,L,M,MULTIP,LAMDA,CO) 
Wire (Gl). . bee 


FORMAT(3X,'a',I1,I1,' = ',F25.12) 
FBAR=FBAR+T**L*CO*DEXP (-LAMDA(I)*T) 
CONTINUE 

CONTINUE 

WRITE(6,15)T,FBAR 

FORMAT (/,1X,'P(T =' 456.2.) =e) 
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40 


50 


60 


70 


90 


STOP 
END 

FUNCTION IFAC(N) 

JFAC='1 

IF(N.NE.O) THEN DO 

DO 40 J=1,N 

JFAC=JFAC*J 

CONTINUE 

END IF 

IFAC=JFAC 

RETURN 

END 

SUBROUTINE COEFF(1I,L,M,MULTIP,LAMDA,CO) 
DOUBLE PRECISION C,CO,LAMDA(10), PRODCT 
INTEGER I,L,M,K,KK,LL,MULTIP(10) 
PRODCT = 1.0D0 

DO 50 J = 1,M 

IF(I.NE.J) THEN DO 
PRODCT=PRODCT*LAMDA (J) ***MULTIP (J) 

END IF 

CONTINUE 

CO = 0.0D0 

Meee Mine (T) - I, 

DO 60 KK = 1,LL 

ieee KK =. 1 

CALL DERIV(I,M,LAMDA,MULTIP,K,C) 

CO = CO + C*PRODCT*LAMDA(I)**(K+L) /IFAC(K) /IFAC(L) 
CONTINUE 

RETURN 

END 

SUBROUTINE DERIV(1I,M,LAMDA,MULTIP,K,C) 
DOUBLE PRECISION A,F,C,B(9),LAMDA(10) 
INTEGER L,J,M,L,IK,MULTIP(10) 

A=.0D0 

F=1.D0 

DO 70 J=1,M 

IF(I.NE.J) THEN DO 

F=F* (LAMDA(J)-LAMDA(I) )***(-MULTIP(J) ) 
A=A+ (MULTIP(J)/(LAMDA(J)-LAMDA(TI))) 
END IF 

CONTINUE 

IF(K.GT.1) THEN DO 

i= 

DO 80 L=1,IK 

B(L)=.0D0 

DO 90 J=1,M 

IF(I.NE.J) THEN DO 
B(L)=B(L)+(((-1)**L)*MULTIP(J)* 


*( (LAMDA(J)-LAMDA (I) )**(-L-1)))*IFAC(L) 


END IF 
CONTINUE 


yee) 


80 CONTINUE 

END IF 

IF(K.EQ.0)THEN DO 

C=F 

ELSE DO 

IF(K.EQ.1)THEN DO 

C=(-1)*A*F 

ELSE DO 

IF(K.EQ.2)THEN DO 

C=((A**2)-B(1) )*F 

ELSE DO 

IF(K.EQ.3)THEN DO 

G= ( (-1)*(A**3 )+3*A*B(1) -B(2) )*F 

ELSE DO 

IF(K.EQ.4)THEN DO 

G= (A**Y -G*A*NRIXB ( lig +4XR AB ( 9) ) + 3% 
*BC1)**2-B(3))*F 

ELSE DO 

IF(K.EQ.5)THEN DO 

C=((-1)*(A**5 )+10* (A**3 )*B(1)-10*(A**2 ) * 
*B(2)-15*A*B(1)=*2+5*A=B(3)+10"B( 1 bigamy 
ELSE DO 

IF(K.EQ.6)THEN DO 

C= (AMG - 15% (M4) *B (1) +20*(A**3)* 

*B(2)- 15*(Ax=2 ) BS) +67AE (4 Ge 
w (A**2 ) ¥ (B ( 1. ) ww” ) -60.DO0*A*B ( 1 ) ve 
*B ( ? ) +15*B ( al ) *B ( 3 ) -15*B ( 1 ) *w*~3+410*"B ( ? ) Weed ) “ER 

ELSE DO 

IF(K.EQ.7)THEN DO 

G= ( ( ~] ) Ww (A**7 ) +2 L KARE KB ( a ) -~35*%ARKGKBR ( 2 ) + 35K AWE ZK 
*B ( 3 ) a 21L*A**2*B(4 ) +7 *AXB ( 5 ) -B ( 6 ) -1O5S*A**3% 
*B ( IL ) WED+ DL ONARNIWB ( 1 ) *B ( ? ) -1O05*A*B ( 1 ) *B ( 3 ) -105* 
*B ( lL ewe 2 eB ( ? ) -7O*A* 
“BRB ( 2 ) WED + LOS*A*B ( lL ) WE Z+DPLEB ( 1 ) *B ( 4 ) +35*B ( 2 ) ¥ 
*“B ( 3 ) )*F 

ELSE DO 

IF(K.EQ.8)THEN DO 

Ge (A**8 -J2E*AKKE XB ( i ) +56 *ARKSKB ( 2 ) -7JO*A**GKB ( 3 ) +56* 
AX*3°B (4) -28%A**2*B(5)+8*A*B(6)-B(7)+210% 

wWAWWEAKR ( 1 ) WED-~-SEQO*KAKRS ZEB ( al ) ¥ 
“B(2)+420*A**2*B(1)*B(3)-168*%A*B(1)*B(4)+28% 

*B ( 1 )*B ( 5 ) = G20*A**2 SB ( 1 ) *3+840*A*B ( al ) WwwIWwR ( ? ) +72?80%* 
KAKKIDRB ( 2 ) wwD-PLO*B ( 1 ) WEDB ( 3 ) -2?80*A*B ( 2 ) *wB ( 3 ) -280*B ( i )* 
*B ( » wD +5 6*B ( 2 )*B ( 4 ) +35*B ( 3 ) #72 4+ 105*B( 1 ) WAG ) “FE 

ELSE DO 

IBC: . 9)THEN DO 

C= ( -A**9+36*A**7*B(1)-84*A7 26 Qop Ee Ile A> 2 Sena) aa 
AREAS "B(4) #842 ‘A%% 3*B (5 )-36*A2*2*B (6 )+9%°A*B(7)-B(8)- 
E37 QRARWORB ( 1 ) WED + 193 Ox AY XB ( 1 ) *B ( ? ) -1760*A**3"B ( 1 ) ¥ 
*B(3)+756"A**2*B(1)*B(4)-252"se (1) BS) Scere eecce 
C=C+1260*A**3*B(1)**3-3780*A**2*B(1)**2% 
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*B(2)-840*A**3*B(2)**2+1890*A* 

*B ( 1 ) WUDKB ( 3 ) +1 260*A**I*NB (2 )*B ( 3 ) -378*B (4 )* 

“B (1) #2 -404*A"B(2)*B(4)+84*B(2)*B(5)-945%A% 

xB ( 1 ) wwOG+D5Q20*A*B ( 1 ) *B ( 2 ) w*D+1960*B ( 1 ) Www ZEB @ ) sss 
“B(3)**2-1260*B(1)*B(2)*B(3)+126*B(3)*B(4)-280*B(2)**3 
C=C*F 

END IF 

ND 

MND LF 

END IF 

END IF 


SENTRY 


Tig, 


section Examples subsection All failure rates are dissimilar 
Shorthand no#mation™ 


EXP(.3)+EXP(.32)+EXP(.4)+EXP(.6)+EXP(.62) 


OutpuE 
TIME = 2.00000 
FAILURE RATE MUTE Tia ee Gey. 

0.30000 i 
0.32000 IL 
0.40000 1 
0.60000 il 
0.62000 il 


POLYNOMIAL COEFFICIENT FOR LAMDA 1 IS 
alO = 248 .000000000000 
POLYNOMIAL COEFFICIENT FOR LAMDA 2 IS 
a20Q = =352,. 1427657142650) 
POLYNOMIAL COEFFICIENT FOR LAMDA 3 IS 
a30 = 101.454545454545 
POLYNOMIAL COEFFICIENT FOR LAMDA 4 IS 
a40 = -70.657 I42 659 ke 
POLYNOMIAL COEFFICIENT FOR LAMDA 5 IS 


a5Q = 54.545454545454 


P(L =e eOoh = 0.998171698511 
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Sipseeei on "All “failure™rates™are identical 


Siewtiana motation 


PegedZ.1) + ... @@EXP(2.1) (10 times) 
Output 
TIME = 520000 
FAILURE RATE MUL eit Cr ry 
22 UO OOS 10 


POLYNOMIAL COEFFICIENT FOR LAMDA 1 ARE 


alO = 1.000000000000 
all = 2. 100000000000 
al2 = 2.205000000000 
al3 = 1.543500000000 
al4 = Oo LU5 47 SV GOC0 
al5 = 0.340341750000 
alo = Ore G2 00 
al/7 = D035 7258657 50 
al8 = 0.009380669484 
al9 = 0.002188822880 
Pa > 3.20) = 0. Sn 77 rren 3017 
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subsection Some sSsimilars.and some®di1ssimagtaeetatlure Goue> 


Shorthand notegeewon : 


+ EXP(6.5)+ ... + EXP(6.5) (5 times) 
+ BRK 4.5)+ .. 2 4 ERP (4s 5) (5 times) 
+ EXP(3.5)+ ... + EXP(3.5) (5 times) 
+ EXP(2.5)+ ... + EXP(2.5) (5 times) 
+ EXP( 1. 5 eee (15D (5 times) 
+ EXP(.50)+ ... + EXP(.50) (5 times) 
Output 
TIME = 8.00000 
FAILURE RATE MULTIPLICITY 

6.50000 5 

4.50000 5 

3.50000 5 

2.50000 5 

1.50000 5 

0.50000 5 


POLYNOMIALS COEFFICIENTS FOR LAMDA 1 ARE 


alO = -0.033655585 55 
all = -0.015984638484 
al2 = -0.002968083834 
al3 = -0, 0002550352777 
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al4 


OF. 


000008633539 


POLYNOMIAL COEFFICIENTS FOR LAMDA 


a20 
a21 
a22 
a23 
a24 


POLYNOMIAL COEFFICIENTS FOR LAMDA 


a30 


a3l 


a32 


a33 
lon 


78824. 
28606. 
4242. 
308 

9 


Sew 2573 
-833405. 
ee DO 2s 

= OZ ete 
Solltscrmle 


407500664440 
965700048750 
082627828356 


.298280647481 
.469913057431 


342698666000 
255639648600 
555234616700 
9271830233870 
852439502142 


Jeo 


Dee 


POLYNOMIAL COEFFICIENTS FOR LAMDA 4 ARE 


a40 
a41 
a42 
a43 
a44 


POLYNOMIAL COEFFICIENTS FOR LAMDA 


a50 
a5l 
a52 
a53 


LODO LOO? . 
Poe) 2 . 
340359. 
-14083. 
4142. 


e/ 360153. 


DLO 0. 


~58447. 
el. 


114320258000 
118591437400 
641443397200 
275902489360 
139971320400 


gee 7 LOO OO 
360064797400 
530039485250 


Lhayei dW) 2 Uejces Sif 


one 


5 ARE 


a54 = -297 8977/7 06625267 


POLYNOMIAL COEFFICIENTS FOR) LAMDASS Ake 


a60 = 4718.810424379766 
a6l = =1562.762326/7477 52 
a62 = 203./49271SGoeere 
a63 = -12.97/7288Zo5ce7 
a64 = 0.350737520646 
PCT > “S23 00)) = — 0'..998940855505 
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